Weak phase stiffness and nature of the quantum critical point in underdoped cuprates 
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We demonstrate that the zero-temperature superconducting phase diagram of underdoped 
cuprates can be quantitatively understood in the strong binding limit, using only the experimental 
spectral function of the "normal" pseudo-gap phase without any free parameter. In the prototypical 
(Lai_ x Sr cc )2Cu04, a kinetics-driven d-wave superconductivity is obtained above the critical dop- 
ing 5 C ~ 5.2%, below which complete loss of superfluidity results from local quantum fluctuation 
involving local p-wave pairs. Near the critical doping, a enormous mass enhancement of the local 
pairs is found responsible for the observed rapid decrease of phase stiffness. Finally, a striking mass 
divergence is predicted at 5 C that dictates the occurrence of the observed quantum critical point 
and the sudden suppression of the Nernst effects in the nearby region. 
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Considering the enormous amount of research activi- 
ties devoted to the problem of high-T c superconductivity, 
it is hardly an exaggeration to regard it as one of today's 
most important unsolved problems in physics. Specifi- 
cally in the underdoped region, it is now commonly ac- 
cepted that the low carrier density in the system neces- 
sarily leads to strong phase fluctuation of the supercon- 
ducting order parameter pQ due to its conjugate nature 
to the number fluctuation. Consequently, the transition 
temperature T c is suppressed significantly below the pair- 
ing energy scale that controls all essential aspects of the 
standard theory of superconductivity [2 . The crucial role 
of phase fluctuation [TJ [3l |4] has recently gained signifi- 
cant support from various experiments [5HZ] in both the 
low-temperature superconducting state and the 'normal 
state' above the transition temperature T c , and is likely 
tied closely to many of the exotic properties [TJ I5HT2] in 
this region. 

Nonetheless, besides this general understanding, sev- 
eral key issues remain puzzling in the underdoped re- 
gion. In spite of an uneventful evolution of the one- 
particle spectral function [11], the superfluid density re- 
duces dramatically near the observed quantum critical 
point (QCP) [13 (at the critical doping S c ~ 5.2% for 
doped La2CuC>4), below which superconductivity ceases 
to exist even at zero temperature. The current consider- 
ation of phase fluctuation pQ would only indicate softer 
phase at lower carrier density, but offer no explanation 
for the complete suppression of superconductivity at zero 
temperature at S < Sqcp- Particularly in La2Cu04, 
5qcp is quite far away from the antiferromagnetic (AF) 
phase boundary, rendering the common consideration of 
competing order unsatisfactory. This vanishing of su- 
perconductivity below 5 CJ the nature of the QCP, the 
dramatic reduction of superfluid density nearby, and the 
controlling factor of the value of 5 C1 all remain challenging 
to our basic understanding. 

Perhaps the most puzzling observation is the sudden 



suppression of the observed Nernst effect at T > T c 
around the same critical doping # C [14]. This indicates 
that not only the long-range phase coherence, but also 
the short-range phase coherence is lost near the QCP, 
a phenomenon unexplainable via simple fluctuation sce- 
nario, for example due to low dimensionality. 

In this letter, we demonstrate that these puzzles can 
be quantitatively understood in the strong pairing limit 
of local pairs of doped holes. We obtain the zero- 
temperature phase diagram with no need for any free pa- 
rameter other than the experimental one-particle spectral 
function of the "normal" pseudo-gap state. A kinetics- 
driven d-wave condensate is found at 6 > 5 C , with a 
largely enhanced bosonic mass, > 40m e . In great 
contrast, ground states consisting of fluctuating p-wave 
pairs are found at S < S c: incapable of sustaining a con- 
densate. At S = 5 C , a mass divergence results from the 
degeneracy of local d- and p-wave symmetry, dictating 
the presence of QCP. Correspondingly, near QCP S < 5 C , 
the diverging mass explains the puzzling dramatic reduc- 
tion of phase stiffness in both long-range and short-range. 
Our study provides a novel yet simple paradigm to the 
behavior of local pairs in underdoped cuprates, and is ex- 
pected to inspire new set of experimental confirmation, 
as well as re-interpretation of existing experimental ob- 
servations. 

Conceptually, a phase-fluctuation dominant supercon- 
ductivity hosts relatively negligible amplitude fluctuation 
of the order parameter at low energy /temperature. This 
implies that the effective low-energy Hamiltonian for the 
charge and pairing channels must have integrated out all 
pair-breaking processes to conserve the amplitude of the 
order parameter, for example, as in the x-y model[15j. 
The higher-energy pairing energy scale should then man- 
ifest itself only through a strong "pair-preserving" con- 
straint of the low-energy Hamiltonian. This is in perfect 
analogy to the replacement of repulsion U of the Hub- 
bard model by a "no double occupancy" constraint in 
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FIG. 1: (Color online) Band dispersion from experiment, our 
fitting (a), and the t-J model[20 (b), and the corresponding 
parameters r, r' and r" (c). (d) Weak doping dependence 
of the effective mass of doped holes in different directions 
indicated by the arrows in the inset. 

its lower energy counterparts, say the t-J model. Con- 
sequently, a new paradigm for the low-energy physics 
emerges at T < T c , which differs completely from the 
amplitude fluctuation of the standard theories. In this 
new physical regime, the detail of the pairing mechanisms 
(AF correlation p~6j [17], spin-fluctuation [18 , or formation 
of bi-polaron[19 ) are no longer essential. Instead, the 
physics is now dominated by the effective kinetic energy 
that controls the phase coherence. Since only one energy 
scale is essential in this regime, the low-energy physics 
should be universal and simple. 

We then proceed to derive the effective motion of the 
local pairs in three steps: 1) obtaining the effective kinet- 
ics of the doped holes from the experimental one-particle 
spectral function in the "normal state" pseudogap phase, 
2) splitting the many-body Hilbert space between un- 
paired and paired holes, and 3) re-represent the effective 
kinetics to reveal the motion and scattering of unpaired 
and paired holes. The resulting pivoting motion of the 
pairs would then lead straightforwardly to all our quan- 
titative conclusions without any free parameter. 

Let's define the irreducible effective kinetic kernel r 
of the low-energy carriers, r = G^ 1 — G _1 (in matrix 
notation), through the measured one-particle propagator 
G and a non-propagating Green's function Gl consisting 
of a single pole at the central energy of the band. The real 
part of r thus controls the propagation of the carriers, 
with offdiagonal elements giving the dispersion, like the 
effective hopping matrix elements. The imaginary part 
of r gives the decay of carriers and will be ignored in this 
study, since only the average motion is of interest. The 
dispersion in the measured G can then be captured by 

H = nuclei* + h.c. (1) 

ii' 

Note that this Hamiltonian is only meant to represent 
the average effective kinetics of the fully renormalized 



one-particle propagator, within the low-energy "pair- 
preserving" constraint. It certainly does not contain in- 
formation of the pairing interaction that connects to the 
high-energy sector, or the scattering process that broad- 
ens the spectral function. Therefore, r is to be distin- 
guished from the "bare" hopping parameter t commonly 
used in the Hubbard or t-J model, as r have fully ab- 
sorbed the effects of interactions and constraints. 

Fig. [TJa) shows the dispersion of the main features 
in the measured spectral functions of (Lai_ x Sr :E )2Cu04 
by angular-resolved photoemission spectroscopy 
(ARPES)[2TJ [22]. Also shown are the dispersion 
obtained by diagonalizing Eq. [I] using r as fitting 
parameters. One notices immediately the well-known 
observation that the dispersion is not rigid against the 
doping level. The close resemblance to the previously 
published t-J model solutions [20] in Fig. [TJb) suggests 
that this strong doping-dependent renormalization orig- 
inates primarily from the competition between the bare 
kinetic energy and the AF interaction. Consequently, 
the resulting fully-renormalized first, second and third 
neighbor kernels, Ti, t' and r" demonstrate strong 
doping dependence as well (Fig. [IJc)). Interestingly, as 
5 decreases, r" is found to increase steadily approaching 
the value of r', and then exceeds r' right at Sqcp- 1 
This is apparently not a coincidence, and reveals an 
important clue to the nature of QCP to be discussed 
below. Since the fully dressed t\ is negligibly small 
at the underdoped regime, understandable from the 
strong antiferromagnetic correlation, it will be dropped 
from our further analysis. As a reference, Fig. [ljd) also 
shows the effective mass of the doped holes, m^, in 
three directions. In agreement with the experimental 
observation [23]. m\ is about constant for 6 > 5.2%. 

Next, we split rigorously the many-body Hilbert space 
into two components [24 : one spanned by the tightly 
bound pairs of nearest neighboring holes, b\- = 4t c jl' 
and the other by the unpaired holes, fj a : 

! 

4 t = fl t + E by ji + [bWbf+...} 

jeNN(i) 

c l = fh- £ by it + [bWbf+...} (2) 

ieNN(j) 

where Yl' denotes "sum only the first non-zero term" 
to avoid double counting [24 . Such a real-space hole 
pair can result from the strong AF correlation [4, 25 , 26 
or the formation of bi-polaron[19 , but is to be distin- 
guished from the real-space singlet pair of electron com- 
monly used in the RVB-like construction|27j. Since our 
study will focus only on the charge and pairing chan- 
nels, for a cleaner presentation, the spin indices and the 
associated sign changes will be suppressed from now on. 
Finally, ignoring the higher order terms owning to the 




FIG. 2: (Color online) Illustration of kinetic processes of a 
pair of holes (filled diamond) to its six allowed destinations 
(open diamonds) under the "pair-preserving" constraint (a), 
through t' (solid lines) and r" (dashed lines). The same is 
alternatively represented by ellipsoids denoting a pair and its 
allowed destinations in the lattice of bond-centered pairs (b). 
The yellow area denotes the region of the 'extended hardcore 
constraint' that excludes other pairs. The dotted square gives 
the convenient unit cell in our calculations. 

low carrier density, and inserting Eq. [2] into Eq. [I] we re- 
represent the carrier kinetics as H = H b + Hf + H b f l24l 
(dropping the redundant fjfj in &Jj/j/J &i'j)> where 

H b = J2Tii>bljbi>j + h.c. (3) 

ii'j 

describes the pivoting motion of the pairs. The remianing 
terms and H b f give the motion of unpaired holes and 
their coupling to the pairs [24 . In essence, this represen- 
tation merely reveals the internal processes embedded in 
the experimentally observed one-particle propagator. 

The resulting pivoting motion of the paired holes is 
illustrated in Fig. [2] Consider a single pair of holes 
(blue and red filled diamonds) located in the fermion lat- 
tice in Fig. [2|a). Under the pair-preserving constraint, 
only three potential destinations (empty diamonds) for 
each hole are allowed. Converting to the lattice of bond- 
centered pairs in Fig. [2^b), one finds two inequivalent 
sites in a square lattice, each of which can hop to four 
first neighboring sites with r', but to only two second 
neighboring sites with r" along their designated direc- 
tion (denoted by the direction of the ellipsoids). As we 
will present elsewhere, this pivoting motion of the pair 
can actually be derived more completely to include both 
the real and imaginary part of r via the equation of mo- 
tion. Alternatively, one can also arrive at the same via 
the ladder diagrams in the strong binding limit in per- 
turbation theory. 

Note that the hole pairs 6's are under a strong 'ex- 
tended hardcore constraint': b]-b],-, = if i = i' or 
j = j' . This is inherited from the Pauli exclusion princi- 
ple of the original fermion operators cj^cj^ = and that 
double occupancy of electrons are not allowed in the low- 
energy sector. In Fig. |5|b), this constraint forbids occu- 
pation of the six potential hopping destinations of each 
pair by another pair. Therefore, these real-space pairs 
are correlated via a "polite" policy that no one is to get 
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FIG. 3: (Color online) The band dispersion of the hole pairs 
without the extended hardcore constraint, at S > ^gcp(a), 
S — ^Qcp(b), and S < Sqcp(c). Insets in (a) and (c) il- 
lustrates the Wannier function corresponding to the lowest 
band, and the kinetic driving force for the local symmetry, 
(d) Strongly enhanced effective mass of the pairs, m*, and 
the doped holes, m^. 

into others' immediate hopping path. This constraint can 
be considered as an infinite short-range repulsion that 
dominates the inter-pair interactions in most cases, and 
is responsible for stabilizing our bosonic system against 
potential phase separation [28 . 

Now, we will show that this rigorously derived, simple 
Hamiltonian gives quantitatively all the above mentioned 
essential aspects of superfluid in underdoped cuprates, 
without any free parameter. For convenience, a unit cell 
is chosen to include four sites as shown in Fig.[2|b). This 
choice explicitly allows one s-, two p-, and one d-wave su- 
perposition within the unit cell, and equates the doping 
level per unit cell in this lattice and that in the standard 
fermion lattice. Fig. [3^a) illustrates the resulting band 
structure of H b at doping level of the superconducting 
phase S = 15% > Sqcp before including the extended 
hardcore constraint. It shows that at low enough tem- 
perature a Bose-Einstein condensate (BEC) would take 
place at a single minimum at momentum q = 0, with a 
pure (i-wave symmetry (red color) . As in standard dilute 
bosonic systems, one thus expects a d-wave superfluid 
with finite stiffness, once given a scattering length de- 
rived from the extended hardcore constraint. 

The (i-wave structure of the pair is better illustrated in 
real space via the corresponding Wannier function, com- 
puted from the Fourier transform of the Bloch functions 
of the lowest band. Fig. [3|a) shows that the low-energy 
pairs has clear d-wave symmetry with nodes along the 
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(tt, 7r) directions of the standard Fermion lattice, in per- 
fect agreement with the experimental observations. [TTJ 
l29l [30] (Notice that our fermionic lattice is rotated by 
45° from the usual convention.) 

We stress that our resulting local d-waxe symme- 
try is completely driven by the fully screened kinetic 
energy. [5U [32] It originates from the dominance of posi- 
tive t' of the local pair, which prefers energetically oppo- 
site sign of the wave function across first neighbors, thus 
favoring a d-wave symmetry (see Fig. [3^a)). In compari- 
son, the positive r" favors opposite sign across the second 
neighbors, thus p- wave symmetry (see Fig.[3jc)). There- 
fore, t' and t" compete by lowering the band energy of 
d- and p-bands, respectively. 

This explains the long-standing puzzle of lack of su- 
perconductivity at lower doping (6 < Sqcp)- Since in 
this region r" > r' [see Fig. [IJc)], Fig. |3jc) shows local 
p-wave pair has lower energy than d-wave pairs. Further- 
more, the structure of the lattice in Fig. |2|b) dictates a 
line of degenerate p-states in momentum space along the 
anti- nodal directions. The pairs can therefore populate 
any arbitrary state along this line without ever forming 
a BEC. That is, the pairs in the system have local p- 
wave symmetry without possible global phase coherence, 
an effect of quantum fluctuation beyond the original con- 
sideration of thermal phase fluctuation pQ. These inco- 
herent p-wave pairs might also offer a new paradigm to 
the insulating glassy behavior of the electronic structure 
in this region, for example the strong scattering in the 
charge channel. 

The competition between d-wave and p-wave also offers 
a natural explanation of the dramatic phase softness and 
the low superfluid density of the underdoped cuprates. 
Indeed, even near the optimal doping (8 « 15%), the 
comparable value of r" and r' leads to a large effective 
mass of the pair m* = (% 2 /l 2 )d 2 €k/dk 2 w 12m^ w 59ra e 
(I being the lattice constant). This gives a rather long 

penetration depth A = \J £^t n ~ 7000A (taking n s ~ 5 
per unit cell), in reasonable agreement with the experi- 
mental value [55] . Furthermore, as S approaches Sqcp, 
t" grows to the value of r', reducing the separation of 
the d-band and the p-band, and in turn flattening the 
d-band. The effective mass of the d band thus increases 
significantly (Fig. [3jd)) . Since the superfluid stiffness 
correlates with the inverse of the effective mass, the large 
mass enhancement dictates the observed very soft phase. 

This analysis reveals the simple yet exotic nature of the 
observed QCP at the end of the underdoped supercon- 
ductivity region Sqcp = 5.2%: It is dictated by the di- 
verged effective mass of the local pairs [Fig.^d)]. At this 
point t' — t" and the d-wave and p-wave pairs become 
locally degenerate and the d-band is thus completely dis- 
perseless, as shown in Fig.|3jb). Since the effective mass 
now diverges [Fig. [3^c)], the pairs can no longer propa- 
gate and allign the phase to develop a condensate even at 
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FIG. 4: (Color online) Experimental supports of predicted 
mass divergence via (a) non-linear doping dependence of in- 
verse penetration depth and (b) non-linear correlation betwen 
inverse penetration depth and transition temperature T c . 

zero temperature. In essence, it is the perfect quantum 
mechanical interference between the first- and second- 
neighbor hopping that renders the local pairs immobile 
and in turn disables the phase coherence of superconduc- 
tivity. 

The diverged effective mass also account for the 
dramatic suppression of diamagnetic response near 
$qcp [34 . Indeed, the diverging effective mass would not 
only suppress the long-range phase coherence but also the 
shorter-range coherence responsible for a strong diamag- 
netic response. The observed dramatic reduction of the 
Nernst signal [14 near the QCP is thus understandable. 

Our predicted mass divergence near QCP is actu- 
ally strongly supported by experimental measurements of 
penetration depth A of the underdoped YE^CusO^ sam- 
ples. Figure [4^ a) shows that over the entire underdoped 
region, the measured A -2 [35] deviates significantly from 
the simple A -2 oc 8 relationship to be expected if the ef- 
fective mass is assumed constant. On the other hand, our 
theory with large doping-dependent effective mass repro- 
duces very nicely the experimental observation. A even 
more direct evidence is provided by the recent measure- 
ment on the extremely underdoped YE^CusO^ sam- 
ples near the QCP [6]. The observed relationship be- 
tween low-temperature A -2 and T c in Fig. [4^b) shows 
a strong non-linear dependence. In fact, the same be- 
havior has also been observed via mutual inductance [36J. 
The zero slope at A — ^ can be interpreted as an in- 
direct evidence of the mass divergence, and our the- 
ory reproduces very nicely the experimental observation. 
(T c is determined by thermally exhausting all the pairs 
J^duDOS(uj)/(e^ T ^ - 1) = 5.) In addition, the mass 
divergence would also make sense the observed dramatic 
enhancement of the residual specific heat [37] and of the 
isotope effect [38 near the QCP. Nonetheless, a direct 
measurement of the superfluid inertia or the effective 
mass of the pair above T c would be highly valuable. 

In conclusion, we demonstrate that all the key feature 
of superconductivity in the underdoped cuprates can be 
described quantitatively in the strong binding limit, with- 
out use of any free parameter. The d-wave symmetry is 
found to originate from the renormalized kinetic energy, 
and the observed superconductivity can be understood 
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as a superfluid of a dilute real-space hole pairs. Our re- 
sult explains the lack of superconductivity at S < 5qcp 
via quantum fluctuation associated with incoherent local 
p-wave pairs. In the underdoped regime, a large effective 
mass enhancement of the hole pairs is found responsi- 
ble for the observed weak phase stiffness. Finally, the 
observed S = 5.2% QCP is found dictated by the diver- 
gence of the effective mass of the hole pairs, which also 
make sense the dramatic reduction of diamagnetic re- 
sponse (the Nernst effect) near the QCP. These successes 
support strongly a simple description of phase fluctuating 
dominant superfluidity for the underdoped cuprates and 
enable further reconciliation of seemingly contradicting 
experimental conclusions in the field. 
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